Mon. Not. R. Astron. Soc. 000, [THZ] (2008) Printed 18 March 2009 (MN WT^ stylo file v2.2) 



Mapping low and high density clouds in astrophysical 
nebulae by imaging forbidden line emission^ 



O : J. E. Steiner^t, R. B. Menezes\ T.V. Ricci^ and A. S. Oliveira^ 



O 



00 



o 



> 

m 

(N 
rn 

rn 
O 

o 



^ Instituto de Astronomia, Geofisica e Ciencias Atmosf ericas , Universidade de Sao Paulo, 05508-900, Sao Paulo, SP, Brasil 
^IP&D, Universidade do Vale do Paraiba, Av. Shishima Hifumi, 2911, CEP 12244-000, Sao Jose dos Campos, SP, Brasil 



Accepted . Received 



ABSTRACT 

Emission line ratios have been essential for determining physical parameters such 
as gas temperature and density in astrophysical gaseous nebulae. With the advent 
of panoramic spectroscopic devices, images of regions with emission lines related to 
these physical parameters can, in principle, also be produced. We show that, with 
observations from modern instruments, it is possible to transform images taken from 
density sensitive forbidden lines into images of emission from high and low-density 
clouds by applying a transformation matrix. In order to achieve this, images of the 
pairs of density sensitive lines as well as the adjacent continuum have to be observed 
and combined. 

Wc have computed the critical densities for a scries of pairs of lines in the infrared, 
optical, ultraviolet and X-rays bands, and calculated the pair line intensity ratios in 
the high and low-density limit using a 4 and 5 level atom approximation. In order 
to illustrate the method we applied it to GMOS-IFU data of two galactic nuclei. 
We conclude that this method provides new information of astrophysical interest, 
especially for mapping low and high-density clouds; for this reason we call it "the 
Id/hd imaging method". 

Key words: Atomic processes - techniques: image processing - techniques: spectro- 
scopic - ISM: clouds. 



1 INTRODUCTION 

Forbidden line intensity ratios from given species (O^, 
O"*"^, N"*", S^) have been widely used in the literature to 
derive average electron temperatures and densities in as- 
trophysical nebulae. The method f or measuring the elec- 
tron t emperature was suggested by iMenzel. AUer &: HebbI 
while the idea of using the [O ll] line intensity 
ratios to measure electron de nsities was suggested by 
[au er. Uffprd fc van Vleckl (|l949l ) and worked out quantita- 
tivelv bv lSeatonI (|l95l ifor both [O ii] and [S ii] li n es. A n 
early review of these methods is given by ISeatonI (|l96d ). 
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(Australia), Ministerio da Ciencia e Tecnologia (Brazil) and SE- 
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Since then, intensity ratios for lines from many other species 
have been proposed and used. 

Electron temperatures are sensitive to the intensity ra- 
tio of the auroral to the nebular components, for example, 
for the Carbon-like 2p^ and the Silicon-like 3p^ ions. A clas- 
sical intensity ratio is that of [Om] /A4363/7A5007. Aver- 
age electron densities are obtained, for example, from the 
Nitrogen-like 2p^ and Phosphorus-like Sp"^ ions. Intensity 
ratios, used very often, are those of [On] /A3726//A3729 
and [S ii] JA6716//A6731 hues. With the development of in- 
frared, ultraviolet and X-ray detectors, other pairs of lines 
have also been used. A comprehensive r eview on the subject 
is given in lOsterbrock fc Ferlandl (|2006l ). 

In the traditional single aperture spectroscopic ap- 
proach, one obtains a single spectrum of a given object and 
only one intensity ratio is measured. This has the obvious 
disadvantage of providing a single average property (tem- 
perature and/or density) for the object along the slit. With 
the development of two-dimensional spectroscopic devices 
such as the Integral Field Units (IFU) and Fabry-Perot in- 
struments, a new approach is possible as one can obtain si- 
multaneously the average property along the line of sight for 
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each point of the object on the sky, if the object is spatially 
resolved. 

In this paper we present a method of transforming two 
images of density-sensitive emission lines into two other im- 
ages, of high and low-density cloud emission. We demon- 
strate the method by applying it to two objects with ex- 
tended nebular emission. We will not discuss the observa- 
tions in great detail nor the physics of each object. The 
reader is referred to the papers in which the objects are dis- 
cussed, for specific information. The focus of this paper is to 
present the method and show some results, for illustration 
only. 



2 EMISSION FROM LOW AND 
HIGH-DENSITY CLOUDS 

The pair of [Sll] AA6716/6731 lines is, perhaps, the most 
used density-sensitive pair of lines. As all forbidden lines, 
these are also sensitive to the gas density. Their ratio 
is density-sensitive because they come from distinct up- 
per levels. As a consequence, the critical density for de- 
excitation of the A;ti(6716 A) line is Nc — 1,400 cm""', 
while for the Ahd(6731 A), is Nc = 3,600 cm~^. The ra- 
tio i?(Ai<i6716/Ahd6731) has a limit of Ru = 1.44 for low- 
density clouds (A^e ~ 10 cm~^) and a limit oi Rhd = OAA for 
high-density clouds (A''e ~ 10'* cm~^). In Table[T]we present 
a list of lines that belong to density-sensitive configurations 
with respective critical densities and limit line ratios. 

The two emission images of the lines, \id, of low criti- 
cal density and \hd, of high critical density, have intensities 
liji^id) and Iij{Xhd)- These two images can be transformed 
in two new images Iij{ld) and Iij{hd) of low and high- 



density cloud emission by the transformation 



Iij{hd) 



Ru 



Rh 



-Rhd 

Rid 



liji^id) 

liji^hd) 



(1) 



Simple algebra shows that this holds by calculating the 
emission properties, using the Rid and R^d from Table [T] It 
is also simple to show that the two transformed images have 
the property that 



liji^hd) = Iij{ld) + lij(hd) 



(2) 



Therefore this transformation can be considered a de- 
composition of the image Iij{Xhd) into the high and low- 
density cloud emission images Iij{ld) and Iij(hd). 

If the observed emission from a given object comes from 
low-density clouds only, then all the emission will be in 
Iij{ld). If we have high-density clouds only, all the emission 
will be in image Iij{hd). However, in practice, we may also 
have clouds with intermediate densities. For these clouds 
part of their emission will be computed in one image and 
part in the other; the proportion will depend on how far the 
actual ratio (and density) is from the two extremes. 

The question that arises here is how to define the "in- 
termediate densities", as this range of densities varies from 
species to species (see Fig.[T]|. We will define a "low electron 
reference density", Nr{ld), as the density for which the line 
intensity ratio decreases the low-density ratio limit by 0.1: 



[ArV]13nm/[ArV]8nm 

[CIII]1907A/[CIII]1909A 

[OIII]88nm/[OIII]52nm 

[SII]6716A/[SII]6731A 

[SiXIII]6.74A/SiXIII]6.69A 




Figure 1. Calculated line ratios as functions of A^e for a sample of 
ions. For optical and ultraviolet wavelengths, an electron temper- 
ature of Te = 10000 K was assumed, while for X-ray wavelength, 
Te = 10^ K. We have only plotted a few ions, for illustration. For 
practical reasons, ions with ratios larger than 2 were excluded. 



In the same way a "high electron reference density", 
Nr{hd), is defined as the density for which the line intensity 
ratio increases the high-density ratio by 0.1: 



R{I{Xid)/Ii\Hd)) = Rhd +0.1 



(4) 



RiIiX,d)/I{Xhd))^ Rid -O.l 



(3) 



We have chosen 0.1 as the tolerance for the line ratios 
somewhat arbitrarily. We could have chosen, for example, 
0.05. The choice has to do with observational uncertainties, 
not with theoretical reasons. Examining Fig. [1] we see that 
these two choices would not produce significant differences 
in the reference densities. Given the typical observational 
errors in lines that could be frequently weak, we believe that 
0.1 is a more reliable quantity than 0.05, for example. 

These reference densities were also computed and are 
listed in Table [T] We will call the low-density regime 
when A''e < Nr{ld), the intermediate-density regime when 
Nr{ld) < Ne < Nr{hd) and the high-density regime when 
Ne > Nr{hd). 

From the data in Table [1] one can see that the low ref- 
erence density is distinct from the low critical density. The 
ratio Nrild) / Nc{ld) ~ 15 — 61 for the infrared, optical and 
ultraviolet lines, but presents a higher range of 39 - 380 
for the X-ray emitting lines. The high reference density is 
similar (within a factor 2) to the high critical density ex- 
cept for the Be-like and He-like ions. For these ions the 
ratio Nc{hd)/Nr{hd) ~ l(f - 10^ and 4 x 10* respec- 
tively. These large numbers are related to the fact that the 
high density lines in these species are intercombination lines 
that have typically much higher critical density when com- 
pared to the corresponding forbidden line. The high refer- 
ence density is, surprisingly, closely related to the low crit- 
ical density for all infrared, optical and ultraviolet lines, as 
Nr{hd)/Nc{ld) ~ 5.1 - 7.3. For the X-ray lines this ratio 
varies from 8.4 to 43. 
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Table 1. Critical and reference densities for selected lines in the infrared, optical, ultraviolet and X-rays, as well as their limit line 
intensity ratios. 



Infrared lines: 



Species 




^hd 




N,{hd) Ria 


Rhd 




Nrild) 


Nrihd) 


Rcfs.'' 




(/im) 


(/im) 


(cm~'^) 


(cm~'^) 






(cm""^) 


(cm"-^) 












C-like 2p^ ions 












Nil 


205 


122 


4.4 xlO 


2.8 xlO^ 1.44 


0.10 




1.8 


2.7 xlO^ 


1,2 


Oiii 


88 


52 


4.9 xlO^ 


3.5 xlO^ 1.78 


0.10 




1.2 xlO 


3.1 xlO^ 


1,2 


NeV 


24 


14 


6.7 xlO^ 


3.5 xlO" 1.10 


0.11 




4.3 xlO^ 


3.5 xlO"* 


1,2 










Si-like 3p^ ions 












Siii 


33 


19 


1.4 xlO^ 


1.2 xlO"* 2.25 


0.08 




2.3 XlO 


9.8 xlO^ 


1,2 


Arv 


13 


8 


2.9 xlO"! 


1.6 xlO^ 1.75 


0.13 




8.3 xl02 


1.9 xlO"' 


1,2 


Optical and ultraviolet lines: 


Species 


^Id 


^hd 


N^ld) 


Nc{hd) Ru 


Rhd 




Nrild) 


Nrihd) 


Refs.'' 




(A) 


(A) 


(cm~'^) 


(cm~'^) 






(cm""^) 


(cm"-^) 












N-like 2p^ ions 












Nl 


5200 


5198 


6.1x10^ 


2.0x10^ 1.49 


0.38 




3.6x10 


3.8x10^ 


1,2,3 


On 


3729 


3726 


9.4x10^ 


4.4x10^ 1.49 


0.26 




4.5x10 


5.5x103 


1,2 


Ne IV 


2424 


2422 


1.8x10"' 


1.3x10^ 1.50 


0.16 




7.6x10^ 


1.2x10^ 


1,2 


Nav 


2068 


2067 


9.5x10"* 


1.1x10*5 1.50 


0.10 




3.8x10^ 


6.9x10^ 


1,3 










P-like 3p'' ions 












Sii 


6716 


6731 


1.4x10^ 


3.6x10^ 1.44 


0.44 




8.1x10 


5.9x10^ 


1,2 


CI III 


5517 


5537 


1.2x10* 


6.7x10"* 1.41 


0.21 




6.0x10^ 


6.6x10"* 


1,3 


Ar IV 


4711 


4740 


7.1x104 


4.5x10^ 1.42 


0.18 




2.9x10^ 


3.6x105 


1,2 


Kv 


4123 


4163 


2.9x10^ 


3.8x10^ 1.37 


0.09 




1.3x10"* 


1.8x10" 


1,3 










Be-like 2s^ ions 












Cm 


1907 


1909 


7.4x10"' 


9.7x10* 1.51 


7.5x10- 


■5 


2.6x10^ 


4.9x10^ 


1,2,3 


NiV 


1483 


1487 


1.4x10^ 


4.4x10^ 1.48 


3.2x10" 


■5 


4.8x10^ 


9.3x10^ 


1,2,3 


Ov 


1214 


1218 


3.7x10^ 


2.6x10*° 1.44 


1.5x10" 


-5 


1.3x10"* 


2.5x10" 


1,2,3 










Mg-like 3s^ ions 










Sim 


1883 


1892 


4.2x10* 


3.2x10*0 1.49 


1.5x10" 


6 


1.5x10^ 


3.1x10^ 


1,2,3 


X-ray lines: 


Species 


^Id 


^hd 




Nc{hd) Ru 


Rhd 




Nr{ld) 


Nr{hd) 


Refs.'' 




(A) 


(A) 


(cm~'^) 


(cm~'^) 






(cm"-^) 


(cm"-^) 












He-like Is-^ ions 












Cv 


41,46 


40,71 


2.2x10^ 


4.4x10*5 12.65 


2.4x10" 


■6 


5.8x10" 


9.5x10*" 


1,4,5 


NVI 


29,53 


29,08 


1.9x10^0 


1.6x10*'^ 5.20 


2.4x10" 


■6 


1.2x10* 


3.3x10** 


1,4,5 


Ovii 


22,10 


21,80 


9.0x10^° 


4.6x10*'^ 3.28 


2.6x10" 


6 


1.0x10^ 


1.0x10*2 


1,4,5,6 


Si XIII 


6,739 


6,686 


1.1x1014 


4.0x10*8 1.91 


3.1x10" 


6 


2.8x10*2 


9.3x10*"* 


1,4,5,6 



References: 1- iRal chenko ct al.l ([20081 ) ; 2- LO stcr brock fc Fcrland l|2006h : 3- lAUer fc CzvzakI l| 19791 ): 4- IPorauet fc Dubaul (|2000l ): 5- 
iPradhan. Norcross fc Hummer ( 198ll): 6- IZhan g & Sampson. Q983). 



3 CALCULATIONS OF LINE INTENSITY 
RATIOS, CRITICAL AND REFERENCE 
DENSITIES 

Pairs of lines other than [S ll] , but from ions with the 
same kind of configuration, include [CI ill], [Ariv] and [K v], 
all in the optical region of the spectrum. These are the 
Phosphorus-like 2p'^ ions. The Nitrogen-like 2p* ions include 
[Nl] and [On] in the optical and [Neiv] and [Nav] in the 
ultraviolet (Table [l]). 

In the infrared, pairs of lines of the Carbon-like 2p2 and 
the Silicon-like 2p^ are also useful density-sensitive lines. 



With the advancement of ultraviolet and X-ray astronomy, 
other types of species have been used, namely, the intensity 
ratio of the forbidden line to the inter-combination lines of 
species such as the Beryllium-like 2s^ and the Helium-like 
Is^ ions. Typical examples are those of the intensity of [C ill] 
/(A1907 A) to that of Cm] /(A1909 A) and [O vii] /(A22.10 
A) to that of O vil] /(A21.80 A), (see Table[l|. 

In order to obtain the line ratios and critical densities 
of the ions, we obtained the populations relative to the up- 
per energy levels involved in the transition. For this we cal- 
culated the equilibrium taking into account the coUisional 
excitation, de-excitation and radiative decay; no recombina- 
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tion was considered. For the C-like, Si-like, N-like and P-like 
ions, we used a five level approximation. For the Be-like and 
Mg-like ions a four level approximation was used. For the 
X-ray He-like ions we used a five level approximation but 
without taking into account the radiative transitions within 
the sub-levels of the triplet. For the Cv and Nvi ions we 
also did not considered the coUisional transitions between 
these sub-levels, as these constants were not to be found in 
the literature. References for the atomic parameters for each 
ion are given in Table [1] We only show ions for which we 
were able to find all atomic parameters necessary to perform 
all the calculations. 

For the infrared, optical and ultraviolet lines we have 
assumed a nebular electronic temperature of T* — 10* K 
while for the X-ray lines we adopted Tx = 10^ K. In real 
situations the temperature may be somewhat different; in 
such cases, for the infrared, optical and ultraviolet lines the 
following equations may be applied 

Nr{T,) = Nr{T*){T,/T*Y (5) 
N,{T,) = NciT*){n/T*f (6) 
and 

R{T,) ^ R{T*){Te/T'y (7) 

Values for these indices are listed in Table [2] for all ions 
for which we found all relevant atomic parameters for at least 
two temperatures. As can be seen in Table [2] the line ra- 
tios Rhd{Te) is temperature independent while RidiTe) also 
is nearly temperature independent, except for the infrared 
lines. Notice that for the X-ray emitting line ratios Rid and 
Rhd, we tabulated the r atios according to the definition of 
iPorauet fc Dubaul (|200(]| '). 

For the X-ray lines, however, the above equations are 
not a good approximation. A better description is given by 

logR{T^) = logR(Tx) + A[{Tx /T,f - 1] (8) 

and 

logNin) = log[N{Tx)] + B logTx/T, + C{logTx/T,f (9) 

Constants for equations [8] and [9] are given in Table [S] 
For the calculations of the constans in this table, we have not 
taken into account the radiative and coUisional transitions 
between the triplet sublevels. 

Where possible, all ratios obtained were double checked 
using the CLOUDY program. One should note that, in our 
calculations, some line ratios are slightly different from those 
provided by the Cloudy models. Examples are the high- 
density limit of the [O ll] line ratio and the low-density limit 
of [Slll] and [Nev] infrared line ratios. For these ratios. 
Cloudy provides a difference of -1-0.10, -1-0.25 and -1-0.15, 
respectively. Our understanding is that these few and small 
differences arise because of recombination contribution to 
these lines and, perhaps, due to different electron tempera- 
tures (see Table [2}. 

4 APPLICATION: THE NUCLEI OF THE 
LINER GALAXY NGC 4736 AND OF THE 
STARBURST/AGN GALAXY NGC 7582 

We have taken two examples of observations from the Gem- 
ini data bank to illustrate the present work. The objects 



Table 2. The o, /3 and 7 indices for the ions for which data are 
available in the literature. 



Ion 


a{ld/hd) 


I3{ld/hd) 


■y{ld/hd) 


Om 


0.52/0.27 


0.34/0.30 


-0.24/-0.03 


Siii 


0.82/0.32 


0.40/0.36 


-0.24/-0.08 


Ni 


-0.50/-0.55 


-0.55/-0.52 


-0.01/0.00 


On 


0.42/0.28 


0.35/0.44 


-0.01/0.00 


Neiv 


0.44/0.42 


0.43/0.45 


-0.03/-0.01 


Sii 


0.54/0.29 


0.44/0.50 


-0.03/0.00 


Cni 


0.27/0.30 


0.31/0.32 


0.02/0.00 



Table 3. The constants for the X-ray line property variations 
with temperature. 



Constant 


Nc{ld/hd) 


Nrild/hd) 


B{Cv) 


-0.198/-0.322 


-0.386/-0.202 


C(Cv) 


0.098/-0.020 


-0.126/0.102 


S(Nvi) 


-0.208/-0.357 


-0.435/-0.185 


C(Nvi) 


0.102/-0.038 


-0.114/0.117 


B(0 VII) 


-0.186/-0.359 


-0.378/-0.164 


C(Ovii) 


0.134/-0.030 


-0.079/0.164 


S(Sixiii) 


-0.147/-0.333 


-0.738/-0.053 


C(Sixiii) 


0.107/-0.012 


-0.223/0.136 


Constant 


Rld/Rhd 




A{Cv) 


0.014/0.027 




<5(Cv) 


1.21/1.00 




A(Nvi) 


0.024/0.033 




5{N VI) 


1.05/1.00 




A(Ovii) 


0.017/0.038 




5(0 VII) 


1.15/1.00 




^(Sixiii) 


0.147/0.072 




<5(Sixiii) 


0.889/1.00 





analysed here are two galactic nuclei suspected for having 
nuclear activity. NGC 4736 is a nearby LINER galaxy while 
NGC 7582 is a starburst galaxy with evidence of being an 
AGN (Active Galactic Nucleus) from its hard X-ray emis- 
sion. They were observed with the purpose of clarifying their 
nature, that is, to detect and characterize a hypothetical 
AGN. In Fig. [2] we display HST images of the regions of 
interest for both galaxies. 

The data were taken with the IFU-GMOS on the Gem- 
ini telescopes. This instrument has 750 fibers in the single 
slit mode (500 fibers on the object and 250 on the sky) and 
operates in the 4000 to 11000 A spectral region. Standard 
procedures for the datacube reduction were used. Before 
applying the transformation matrix (equation [T]), we have 
de-convolved the datacube in the spatial dimension with 
a Richardson-Lucy algorithm, using a gaussian PSF with 
FWHM equal to the observed seeing and performing 6 iter- 
ations. The goal of this paper is to show that the proposed 
transformation of images is a useful tool for astrophysical re- 
search; we do not intend to discuss the details of the objects 
themselves. 

In NGC 4736, the extracted images of [S 11] emission 
were achieved with an artificial filter matching the individual 
line width. An adjacent stellar continuum was subtracted. 
No velocity information was attempted, as the emission is 
weak. 
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Figure 2. HST observations of the nuclear region of NGC 4736 {left, in a UV band, F220W, centred at 2300 A) and NGC 7582 {right, 
in a NIR band, F160W, centred at 16,060 A). The scale and orientation of the fields are the same as observed with the Gemini Telescope 
(Figs. |3] and HJl. 



The images show condensations of low (Fig. [3};) and 
high (Fig. [3}1) density gas emission. These two maps form 
the NLR (Narrow Line Region) of the AGN. In this case, the 
location of the BLR (Br oad Line Region) is kn own to be at 
X — +0.2 and y — +0.2 jSteiner et ahllin pressh . Two blobs 
of high-density clouds, separated by 0.6 arcsec are seen near 
the BLR, while only one is seen in the low-density cloud 
map. A second, but much weaker blob is seen 2.4 arcsec to 
the South. 

Interestingly, the double structure that seems to emerge 
in the hd image resembles the doubl e radio source observed 
bv lKording. Colbert fc Falckd (|2005l ). The radio sources are 
separated by about 0.9 arcsec and have an alignment with a 
position angle of about —43°. This can not be associated to 
the two blobs in hd image as they present a position angle 
of PA = -10°. 

In the datacube from the nucleus of NGC 7582 (Fig.H]), 
the [S ll] line emission is significantly stronger. Therefore we 
extracted the images with a similar procedure as above, but 
with narrower filters (co-adding images of contiguous wave- 
length pixels), so that velocity information is also registered. 
Three velocity bands were defined. Red is associated with a 
filter from to +80 km s~^; green is a filter from to —80 
km s^^ and blue is from —80 to —243 km s~^. The blue- 
shifted band is necessary because of the blue wing seen in 
most of the forbidden lines. 

The imag es display, first of a ll, a clear (known) galactic 
rotation map jMorris et al.lll985h . In addition, in the high- 
density image (Fig. |311) the bulk of emission comes from two 
Hll regions. These two Hll regions are associated to the 
starburst character of this galaxy and their emission lines 
overshine that of the AGN. 

The low-density emission image has a complex and 
spatially extended structure. This complexity is due to 
knots that correspond to H II regions and a few hot 
spots, presumably Supernovae Remnants - SNR; in addi- 



tion, a diffuse bluish emission that comes from a wind, pre- 
sumably photoionized by the AGN bu t not as much ob- 
scured, delineating the ionization c one (|Morris et al.lll985l : 
IStorchi-Bergmann fc Bonattdll9^ ) . A paper with a full dis- 
cussion of the datacube of this object is in preparation. 



5 CONCLUSIONS 

The density sensitive emission maps should not be confused 
with maps of line ratios I(A6716/A6731), which are also used 
in the literature. Such an image does not indicate the in- 
tensity and is frequently very noisy when the denominator 
becomes small or even zero. 

The reader should be cautioned to the fact that if a 
given spatial pixel coincides with a low-density cloud and a 
high-density cloud, projected along the same line of sight, 
one may be confused as this could be interpreted as emission 
coming from an intermediate-density cloud. If velocity sepa- 
ration is possible, as we show in the example above, then the 
chance for the two clouds to coincide in spatial projection 
and simultaneously in velocity, diminishes. 

Although the sum of the images of high-density and low- 
density is mathematically equal to the image of the line of 
highest critical density (see equation [2| , this is not always 
obvious in the displayed images because of distinct LUT 
(look-up table) dynamical ranges and choices. 

For the present methodology to be applicable, the cloud 
velocity distribution must be smaller than half the separa- 
tion of the two lines, otherwise there will be superposition 
of emission on the same spectral pixel. This is usually the 
case for planetary nebular, Hll regions, symbiotic stars etc., 
but may not be the case for novae and supernova remnants 
and most of the active galactic nuclei. It is worth mentioning 
that the various ionic configurations have very distinct line 
separations. For instance, the infrared lines are quite sepa- 
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Figure 3. The nuclear region of the LINER galaxy NGC 4736: (a) Top left: observed imago of the [Sll] emission Iij{Xid 6716 A); (b) 
Top right: observed image of the [Sll] emission lijiXhd 6731 A); (c) Bottom left: transformed imago of low-donsity cloud emission; (d) 
Bottom right: transformed imago of high-donsity cloud omission. 



rated (in velocity) as are the X-ray lines (whose separation is 
~ 5000 km s^^) where this warning does not apply. But for 
the nitrogen-like 2p^ configurations, the line separations are 
typically only 2 A. This is the case for [On] A3736/A3729. 
In such situations the technique applies only if the involved 
velocities are smaller than ~ 100 km s^^. 

Two-dimensional optical devices such as Integral Field 
Units (IFUs) and Fabry-Perot have advanced significantly in 
recent years; the method proposed in this paper, when com- 
bined with such instruments, may provide a powerful tool 
for science. Future developments may extend this method 
to other wavebands, so that many of the pairs of lined listed 
in Table [1] could become potentially useful. 

We conclude that this "Id/hd imaging method" may be 
useful in mapping low and high-density cloud distribution 
in astrophysical nebulae. In addition to density maps it also 



provides a density-related intensity image in specific forbid- 
den lines. 
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Figure 4. The nuclear region of the AGN/starburst galaxy NGC 7582: (a) Top left: observed image of the [Sll] emission Iij{Xid 6716 A) 
(b) Top right: observed image of the [S ll] emission Iij{Xhd 6731 A); (c) Bottom left: transformed image of low-density cloud emission 
(d) Bottom right: transformed image of high-density cloud emission. The RGB colors have the following ranges: R: from to 80 km s~^ 
G: from to —80 km s~^ and B: from —80 to —243 km s~^. Image (d) is dominated by two Hll regions centred at (x = 0.7; y = 0.5) and 
(—0.2; —0.5) while image (c) is dominated by the H II regions, the ionization cone and a few hot spots of unidentified origin, perhaps 
SNRs. 
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